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We report results of a Monte Carlo simulation of the i/i* quantum field theory using multigrid simulation 
techniques and a refined discretization scheme. The resulting accuracy of our data allows for a significant test of 
an analytical approximation based on a variational ansatz. While the variational approximation is well reproduced 
for a large range of parameters we find significant deviations for low temperatures and large couplings. 



1. INTRODUCTION 

Among the methods to treat one-dimensional 
quantum systems analytically the variational ap- 
proximation has been shown to be a very useful 
one [1-3]. Since, however, it is hard to give precise 
inherent error estimates it is desirable to check 
the variational method against independently ob- 
tained numerical data. The difficulty is to achieve 
sufficient accuracy. Standard path-integral simu- 
lations suffer from well-known draw backs, nota- 
bly systematic errors due to the necessary dis- 
cretization and severe slowing down in the conti- 
nuum limit. 

In view of these difficulties it is therefore gra- 
tifying that recently some algorithmic improve- 
ments developed for spin systems and lattice field 
theories could successfully be transferred to path- 
integral simulations [4]. Multigrid simulation 
techniques in particular have been shown to elimi- 
nate slowing down in the continuum limit for one- 
particle systems [5]. Here we will present Monte 
Carlo (MC) simulation data for the (j>^ quantum 
chain obtained by combining multigrid techniques 
with a refined discretization scheme. Finally we 
take the spatial continuum limit of this system. 

2. THE MODEL 

The quantum field theory in its spatially 
discretized version is a system of N harmonically 
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coupled anharmonic oscillators described by the 
partition function 
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where the potential is given by 
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Here /3 = l/ksT denotes the inverse tempera- 
ture, (fii = d4>ildu, and X'[{(/)j(m)}] is the usual 
path-integral measure. We always employed pe- 
riodic boundary conditions, (/)o = 4'N- Following 
the notation of Ref.[3] we introduce a dimension- 
less coupling constant Q = TiLOi/ Eg and parame- 
ter R = LOo/uJi, where Eg = (2/3)j4awo'^i is the 
energy of the classical static kink and R its length 
in units of the lattice spacing a. The reduced tem- 
perature will be denoted by t = ksT/Eg. 

The variational approach [1-3] starts from a 
quadratic trial Hamiltonian with free parameters 
which are determined by optimizing the Jensen- 
Peierls inequality for the free energy. Already for 
one-dimensional quantum systems [1,3], however, 
the resulting set of coupled equations is extremely 
complicated and explicit solutions have been ob- 
tained in the literature only for the limiting ca- 
ses of high and low temperatures and for small 
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coupling Q. As a result one obtains an effec- 
tive classical partition function which can be eva- 
luated, e.g., by standard transfer integral techni- 
ques. The small coupling limit, given in eq.(3.6) 
of Ref.[3], seems to be the most useful one and 
will be used for comparison with our MC data 
for the internal energy u per site and the specific 
heat c per site. More precisely, we will be inte- 
rested only in the anharmonic contribution Rdu 
resp. Rdc to these quantities. 



3. SIMULATION TECHNIQUES 

For the numerical work the partition function 
(l)-(3) was further discretized in euclidean time 
direction using the Takahashi-Imada scheme [6] . 
This guarantees a convergence of order e^*, where 
e = [3 1 L and L is the Trotter number. For the 
evaluation of the anharmonic contributions the 
discretization error was further reduced by sub- 
tracting the exact values for the harmonic con- 
tribution at finite Trotter number L which can 
readily be derived also for the Takahashi-Imada 
scheme [8]. 

Since for local update algorithms we expect a 
quadratic slowing down in the continuum limit of 
small e [5] we applied at each site a multigrid W- 
cycle with piecewise constant interpolation along 
the Trotter direction using single-hit Metropolis 
updating with one presweep and no postsweep. 
The Metropolis acceptance rates were adjusted 
to be 40 — 60% on the finest grid and the same 
step width was used for all multigrid levels. 

For various Q and t we measured the internal 
energy using an optimal combination of the "ki- 
netic" and "virial" energy estimators [7,8] with 
measurements taken every second sweep in runs 
of 400 000 sweeps after discarding 2 000 sweeps for 
thermalization. The specific heat was measured 
by simple numerical differentiation of our com- 
bined estimator. All statistical errors were com- 
puted by jackkniving the data on the basis of 500 
blocks. The energy measurements are practically 
uncorrelated since the integrated autocorrelation 
time for both estimators was always bounded by 
r'"' < 4. Within this bound r'"' tended to be 
larger for low t and large Q. 
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Figure 1. Anharmonic contributions to the inter- 
nal energy. Solid (dashed) lines show the varia- 
tional approximation for N = 300 {N = cxd). 



4. RESULTS 

Figure 1 shows the measured anharmonic con- 
tributions to the internal energy. Here the num- 
ber of oscillators was N = 300 except for t = 0.05, 
0.30, 0.35, and 0.40 where we simulated a chain 
of # = 200 oscillators. The Trotter number was 
L = 16 for t > 0.20 and L = 32,64,128 for 
t = 0.15,0.10,0.05. By performing a number of 
simulations with smaller Trotter numbers we con- 
vinced ourselves that these values suffice to re- 
duce the remaining discretization to be at most 
of the order of the statistical error. 

Comparing our data with the variational ana- 
lysis we find that the approximation is fully con- 
firmed for high temperatures t and small coup- 
lings Q. For lower t we still find a satisfactory 
agreement if we also take into account finite-size 
corrections. The situation is different, however, 
for low temperatures and large couplings where 
we observe significant deviations from the varia- 
tional approximation. Note that the error bars for 
the data in Fig. 1 are smaller than the data sym- 
bols. For the worst case, Q = 0.6 and t = 0.05, 
the variational approximation deviates from the 
MC results by 56 statistical error bars, even if 
finite-size corrections are fully taken into account. 

Figure 2 shows the measured anharmonic con- 
tributions to the specific heat per site. Since 
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Figure 2. Anharmonic contributions to the spe- 
cific lieat per site. Solid (daslied) lines show the 
variational approximation for N = 300 (N = cxd). 
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Figure 3. Anharmonic contributions to the spe- 
cific heat per site of Fig. 2 rescaled to the conti- 
nuum field theory. 



the estimation of the specific heat involves a dif- 
ference of statistically fiuctuating variables the 
resulting statistical accuracy is greatly reduced 
compared to the estimation of energies. There- 
fore our data for the specific heat do not allow 
for a significant falsifying test of the variational 
approximation, in particular for large couplings 
and low temperatures. 

Let us finally consider the renormalized conti- 
nuum theory. Following the procedure in Ref.[9] 
we have rescaled in Fig. 3 our data in order to 
be able to compare them with the variational ap- 
proximation for the continuum case, i.e., the (j>^ 
quantum field theory in 1 -|- 1 dimensions. For 
small couplings we find again good agreement 
between the variational approach and our simu- 
lation data. 

5. DISCUSSION 

Employing refined path integral MC techniques 
to the quantum chain we find very good agree- 
ment with a variational approximation for small 
couplings and high temperatures. For large coup- 
lings and low temperatures, on the other hand, we 
observe clear discrepancies which can be caused 
either by the small coupling expansion used to 
approximate the effective classical partition func- 
tion or by an inherent failure of the variational 



approximation itself (at this level of accuracy). 
It would therefore be interesting to investigate if 
higher-order corrections in the coupling parame- 
ter Q or higher-order corrections to the variatio- 
nal approach itself can account for the remaining 
discrepancies. 
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